Depletion effects and loop formation in self-avoiding polymers 
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Langevin dynamics is employed to study the looping kinetics of self-avoiding polymers both in 
ideal and crowded solutions. A rich kinetics results from the competition of two crowding-induced 
effects: the depletion attraction and the enhanced viscous friction. For short chains, the enhanced 
friction slows down looping, while, for longer chains, the depletion attraction renders it more frequent 
and persistent. We discuss the possible relevance of the findings for chromatin looping in living cells. 
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The kinetics and thermodynamics of the folding of a 
flexible polymeric chain into a loop are central issues in 
polymer physics 0, S IS 0- Renewed interest in this 
classic problem has been fuelled by the introduction of 
novel manipulation techniques [f| that provide unprece- 
dented insight into the mechanics and flexibility of vari- 
ous biopolymers. In particular, it has been shown that, in 
the cell nucleus, DNA regions separated by several /im's 
on the genetic map can nevertheless be in molecular con- 
tact . Is such looping generated by active mechanisms 
(e.g. molecular motors), or merely by passive thermo- 
dynamic mechanisms (e.g., diffusion)? Some light can 
be shed on these issues by comparing experimental ob- 
servations with the statistics and dynamics of looping 
predicted by general polymer models. The systematic 
application of this strategy has so far been hindered by 
the dependence of cyclization dynamics on many time 
scales even for the simplest phantom polymer models @ . 

Here, we go beyond the treatment of phantom chains 
and focus on the impact of steric effects. We not only 
consider the polymer self avoidance 0, 0] but also incor- 
porate excluded volume interactions of the chain with a 
surrounding crowded environment, treated as a collection 
of small monodisperse globular particles (microspheres) 
which induce an entropic attraction on larger objects in 
solution s]. To date, the investigation of this intrigu- 
ing depletion effect, has been mainly studied in poly- 
dispersed colloidal solutions 0, Q . Understanding how 
crowding affects the behaviour of a single self-avoiding 
polymer thus represents a novel and important topic in 
macromolecular physics. It also has immediate implica- 
tions in cell and systems biology, as cells are so crowded 
with globular proteins and RNAs [H IH IH ■ We 
show that crowding affects the occurrence and persis- 
tence of loops in self-avoiding polymers in diverse ways, 
according to the length of the polymer chain and the 
size of the constitutive monomers. Besides uncovering 
new physics, our results may be relevant to the under- 
standing of chromatin looping in vivo. Specifically, we 
will discuss whether and to what extent the depletion at- 
traction may explain existing observations in cell biology 
that active polymerases, attached to chromatin, cluster 



into supramolccular "factories" of up to fim size during 
transcription and replication |l3j . 

A traditional string-and-bcads model will be used to 
describe polymer chains. We consider two simple cases: 
a plain self-avoiding string of N equally-sized spherical 
beads of radii Jii = i ) ... ) iv=12.5nm, and one with larger 
beads at the ends. This sphere size was chosen to re- 
late our polymer to an eukaryotic chromatin fiber, whose 
effective diameter and persistence length are both ~ 25 
nm 0| . The case of larger end beads is instead moti- 
vated by the study of chromatin loops with an attached 
genome-grabbing machinery (a large transcription or a 
replication complex which locally increases the effective 
fiber diameter [l3j). The potential energy, when the cen- 
tres of beads i and j are at a distance d% j , is 
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where ei and £2 are respectively 0.24 and 70 units of 
thermal energy, kbT, a — 4 nm" 1 , and d® j = Ri + Rj 
is the contact distance of beads i. The first term in eqn. 

enforces the hard-core repulsion for contacting pairs, 
while the second provides an attraction between consecu- 
tive beads on the chain. Interplay between the two terms 
produces a self-avoiding FENE chain where, at tem- 
perature T = 300 K, the distance between consecutive 
beads fluctuates by about 0.5 nm around 25 nm. The 
microspheres have radius r = 2.5 nm and occupy a frac- 
tion (j> = 0.15 of the total available volume. These values 
conservatively reflect the crowding of the cellular envi- 
ronment mostly due to RNA and proteins [lflj. Because 
the value of <j> considered here is moderate (see Fig^i). 
we can resort to the approximate Asakura-Oosawa (AO) 
treatment which does not require to simulating ex- 
plicitly the dynamics of microspheres. More precisely, in 
addition to the term of eqn. Q, the polymer is subject 
to the following effective interaction potential: 
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(c) Average crossing number 



FIG. 1: [Color online] (a) Sketch of a chain of N = 8 beads 
(end spheres highlighted) of radius R — 12.5nm. Micro- 
spheres (r = 2.5nm and A = 0.15) within 7.5 nm of the chain 
surface are shown, (b) Time evolution of the end-to-end dis- 
tance, di t N, in chains of N = 5 beads subject to the effective 
depletion potential accounting for the presence (A = 0.15) 
and absence (A = 0) of the microspheres. Looping occurs 
when di t N < 2(R + r) = 30 nm (shaded area). Probability 
distribution of (c) average crossing number and (d) chirality 
for A = and A = 0.15. 



step function O ensures that the AO depletion interaction 
vanishes at distances > d Q + 2r Q. As r is sufficiently 
smaller than the radii of the chain beads it is legitimate 
to disregard in three- and many-body interactions. 
The evolution of the system, carried out for various TV, 
sizes of end spheres and values of (j) was described by 
overdamped Langcvin dynamics 



li ±f = -d(v c + v AO )/dxf +e?(t) 
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where a runs over the Cartesian components, Si is the po- 
sition of the ith bead and the stochastic white noise term 
obeys the fluctuation-dissipation condition: (£) = 0, 
(£f(*)<f (*0) = 28a,f>8i,jS(t,t) k B T lt . The friction term 
7i was obtained from the Stokes-Einstein Q relation- 
ship: 7l = 67177(1 + 2.5</>)i? 4 where r) = 5 cPHjJ. The 
Langevin equation was integrated numerically by means 
of a predictor-corrector scheme jlfij and a time step of 15 
ps. The viability of eqn. (JSJ) was ascertained by a prelim- 
inary successful comparison of various dynamic and equi- 
librium properties with those produced by underdamped 
dynamics (with masses deduced from typical densities of 
biopolymers, p = 1.35 g/cm 3 (l7jl. 

The dynamical evolution was followed starting from 
randomized non self-intersecting configurations of chains 
with 3 < N < 30 beads. When investigating how looping 
dynamics is affected by the size of the contacting spheres, 



we set N = 10 and varied the radius of the end spheres, 
Ri = Rn, within 12.5 and 43.75 nm. The formation 
of loops was detected by monitoring the end-to-end dis- 
tance, d\^M and comparing it to the range of the depletion 
attraction, di,jv < [d\ N + 2r). To have a well-defined 
comparison term, the same criterion for loop formation 
was adopted in the absence of crowding effects/agents 
(i.e. 4> = 0). Typical evolutions of the end-to-end dis- 
tance are illustrated in Fig. {lfa. The trajectories were 
analysed to highlight how depletion interactions affect 
the looping kinetics and the chains equilibrium struc- 
tural properties. For the latter issue several geometri- 
cal descriptors were considered: the radius of gyration, 
virial coefficients and the distribution of local chiralities 
Xi = w t +2.i+3 • A u i+ i. l+2 ), Ui.i+i being the nor- 

malised bond vector joining residues i and i + 1. For 
looped configurations we also calculated the writhe and 
crossing number, averaged over hundreds of randomly 
oriented two-dimensional projections |l8j . Concerning 
the looping dynamics we instead characterize the evolu- 
tion of the system in terms of the mean first looping time 
(MFLT) and the mean-first unlooping time (MFUT), 
which we establish through the following novel proce- 
dure apt for numerical implementation. For very long, 
and hence thermalised, trajectories the MFLT is obtain- 
able by picking at random unlooped conformations and 
measuring the time to the first looping event. Configu- 
rations in a specific "unlooped time interval" of duration 
t u will be picked with weight equal to t u and their av- 
erage first looping time will be t u /2. The MFLT can 
thus be expressed in terms of the average duration of un- 

looping intervals and its second moment: MFLT= ^-^y- 
An analogous formula holds for the MFUT . The aver- 
age values of MFUT and MFLT (we have verified that 
the first two moments of looping and unlooping intervals 
are finite) and their uncertainties were calculated over 
10 independent trajectories having maximum duration 
ranging from 0.1 s for N — 3 to 10 s for N = 30. 

We first discuss the structural differences of the gener- 
ated configurations. The short-range depletion attraction 
produces a reduction of the effective size of the polymer. 
For the largest N considered, where most conformations 
are unlooped, the radius of gyration decreases by 10% 
when <j) goes from to 0.15. The depletion attraction also 
impacts on the chain structural organization. Indeed, the 
geometrical complexity of looped chains is enhanced by 
crowding and the average crossing number increases with 
4> (Fig- IDO- This is akin to what occurs in random rings 
upon compactification |l8j . At variance with this case, 
the development of a striking trimodal character of the 
chirality distribution indicates the emergence of a pecu- 
liar structural organization (see Fig. QJl). Though the 
chiral biases are local, and hence do not lead to long- 
range structural order, they provide qualitative support 
to the recent suggestion of Snir and Kamien that dcple- 
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FIG. 2: Mean first unlooping time as a function of (a) chain 
length and (b) radius of terminal spheres. The curves in (b) 
correspond to exponential and linear fits of the data, (c) 
Mean first looping time as a function of chain length (points 
and dashed line). The asymptotic law for relaxation time 
in Rouse chains, tr cx N 2u+1 , is shown for comparison, (d) 
Ratio between MFLT in the presence (0 = 0.15) and absence 
(4> = 0) of crowding agents. The typical error is 5%. 



tion effects may be sufficient to drive the formation of 
optimal helices in thick biopolymers 0, . 

We next turn to the analysis of the looping kinetics 
and discuss the behaviour of the MFUT. In the presence 
of microspheres, the two contacting ends arc subject to 
the depletion attraction and it may be anticipated that 
the MFUT is larger than for = 0. This expectation, 
qualitatively perceivable from Fig. \T]p, is confirmed and 
quantified in Fig. which portrays a parallel trend of 
unlooping time as a function of N for = and = 0.15. 
Moderate crowding is enough to increase the time the 
ends spend together (forming a loop) by a factor of 3. 

The behaviour of the MFLT (Fig. Eb and d) is more 
intriguing and harder to anticipate, owing to two oppos- 
ing tendencies. On one hand, as we saw earlier the de- 
piction attraction reduces the effective polymer size and 
hence favours looping. On the other, crowding augments 
the effective viscosity of the medium, thereby slowing the 
diffusive encounter of terminal beads. More precisely, the 
Stokes-Einstein formula gives a 38 % increase of friction 
coefficient when = 0.15 compared to = 0. The re- 
sulting balance between the two opposing effects can be 
established by considering the asymptotic expression for 
relaxation times in Rouse chains with excluded volume. 
The slowest relaxation time in such chains (assimilated to 
the MFLT HQ) increase as jb 2 N 2 " +1 , where v « 0.6 
is the scaling exponent for self-avoiding polymers and 
b is the effective size of the chain monomers estimated 
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FIG. 3: Theoretical behaviour of the ratio between MFLT at 
finite volume crowding, 0, and the MFLT (0 = 0). 



by calculating the second virial coefficient accounting for 
the depletion attraction of of cqn. [3 1]. Indeed, the 
data collected within the explored range of N, appear 
well compatible with this asymptotic relationship (sec 
Fig- lit). The Rouse scaling formula can hence be used 
to quantify how the crowding-induced changes in 7 and b 
ultimately affect the MFLT for large N. For the specific 
case considered here, = 0.15, one finds that the MFLT 
is decreased by about 17% compared to the = case. 
This crude asymptotic estimate is in fair agreement with 
the simulation results for N =15-30 which indicate that, 
when = 0.15, crowding decreases the MFLT by 10 % 
or more. It is interesting to consider, within the pre- 
vious approximate analysis, how the MFLT depends on 
0. This information, obtained from the known functional 
dependence on of the friction and second virial coeffi- 
cients, is shown in Fig. [3] which portrays an intriguing 
non-monotonic dependence: for small crowding < 0.1, 
the diffusive slowing dominates, while larger facilitate 
looping via depletion-induced crumpling of the chain. 

Wc finally discuss the dependence of looping on the size 
of the terminal beads. We considered chains of N = 10 
beads where R\ = Rn were varied between 12.5 and 
43.75 nm. We found that the increase in end spheres 
size had minor effect on the MFLT which, for = 0.15, 
changed by less than 10 % over the explored range of ter- 
minal radii. The near constancy of the MFLT is notewor- 
thy since the change of R\ = Rn from 12.5 to 43.75 nm 
implies a two-fold increase of probability that terminal 
spheres contact internal ones. In contrast, the MFUT 
is greatly affected by changes in terminal radii. As il- 
lustrated in Fig. |2b for = 0.15 it increases approxi- 
mately exponentially as a function of R\ and Rn- This 
behaviour would be expected if the terminal bead exited 
the depletion well (whose depth increases linearly with 
i?i = Rn) according to simple Arrhenius kinetics. How- 
ever, the increase in MFUT also reflects the decreased 
diffusivity of terminal beads as a result of the linear in- 
crease in R of the friction coefficient, 7. This second 
effect, dominated by the former in the presence of crowd- 
ing, is readily visible in the curve pertaining to = in 
Fig. [2). The dotted line represents a fit to the 6 data 
points with a linear relationship in R± = Rn (with rela- 
tive x 2 equal to 1.5), which would be appropriate if the 
diffusion coefficient of the terminal spheres were the only 
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factor slowing unlooping. By carrying out simulations 
for selected values of <j> it was found that increasing <j> en- 
hances the magnitude of these effects. For example, with 
cf> = 0.3, which may also be appropriate for crowding in 
nucleoplasm [l^, the MFLT is ~ 3 times smaller, and 
the MFUT (with R x = R N = 25 nm) ~ 9 times longer 
than with <fi = 0.15. 

We now discuss the possible biological implications of 
our results. As mentioned previously, a large number of 
active DNA (and RNA) polymerases can attach to spe- 
cific segments of a chromatin fiber thus increasing its 
local thickness. Experimental observations have shown 
that these scattered groups of polymerases eventually 
cluster into replication (or transcription) factories thus 
looping the intervening genome. In the case of transcrip- 
tion factories in cukaryotes, these structures contain ~ 10 
or more polymerases and having size ranging from ~ 100 
nm up to ~ fj,m I s there a simple physical mech- 

anism leading to their establishment? Our study sug- 
gests that the formation and persistence of these loops 
can be aided by cellular crowding. First, using a conser- 
vative estimate of = 0.15, we find that the depletion 
self-interaction of the fiber thermodynamically facilitates 
looping (see e.g. Fig. QJ. For example, analysis of end- 
to-end distances in equilibrium shows that bridging the 
two ends of a 750-nm fiber, which would contain 75 
kilobases of DNA 01, costs more than 8 i^T in the ab- 
sence of crowding agents, but less than 7 fcgT in their 
presence. We find crowding diminishes the looping cost 
by 1-2 ksT for all lengths simulated. All, or most, of 
this looping cost may be overcome by the extra deple- 
tion attraction between the thicker ends of the loop, to 
which the transcription machinery is attached, consis- 
tently with recent theoretical predictions ^jj. Second, 
the results of Fig. |21i demonstrate that crowding can also 
aid looping dynamically by facilitating the diffusive en- 
counter of the ends. This intriguing result is supported 
by theoretical arguments summarised in Fig. |3 which 
suggest the effect is robust. Furthermore, crowding sta- 
bilizes loops once the two ends have met. The fact that 
the MFUT has an approximately exponential dependence 
on the height of the depletion well (Fig. [Sja) underscores 
the role that crowding has for the formation of replica- 
tion/transcription factories. Via eqn. we estimate 
the unlooping time of e.g. two transcription complexes 
of radius R w 40 nm with <fr ra 0.3 can easily exceed 
0.1 s - a macroscopic time-scale, ~ the experimentally 
measured persistence time of factories. In living cells the 
life-time of such aggregates is probably longer because, 
besides other physical-chemical factors, many complexes 
may come together (rather than the two considered here) 
and interactions of one large bead with two or more oth- 
ers cooperatively increases the stabilization due to the 
depletion attraction [l3| . 

In conclusion, we have considered various kinetic and 
thermodynamic aspects of polymer looping in a crowded 



medium. The process of loop formation is controlled by 
two opposing effects. On one hand looping is entropi- 
cally aided in crowded media by the depletion effect. On 
the other, the enhanced friction of the medium hinders 
the diffusive encounter of the chain ends. The balance 
of the two effects is found to depend both on the length 
of the polymer chain and on the size of the contacting 
ends. Specific model parameters have been used to show 
quantitatively that crowding-enhanced looping forma- 
tion/persistence may be actually exploited in living cells 
to promote the contact of actively replicated/transcribed 
chromatin regions as observed in recent experiments . 
The approach outlined here demonstrates the viability of 
computational and analytical approaches to investigate 
the novel and stimulating problem of crowding effects on 
a single polymer. It would be interesting to confront the- 
ory and experiments on looping and unlooping times ob- 
tained from single molecules experiments with and with- 
out crowding agents. 
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